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Abstract 

By analyzing the dielectric non-linearity with the Landau thermodynamic expansion, we find 
a simple and direct way to assess the importance of the eighth order term. Following this ap- 
proach, it is demonstrated that the eighth order term is essential for the adequate description of 
the para/ferroelectric phase transition of BaTiOs. The temperature dependence of the quartic 
coefficient /3 is accordingly reconsidered and is strongly evidenced by the change of its sign above 
165°C. All these findings attest to the anomalously strong polarization anharmonicity of this 
material, which is unexpected for classical displacive ferroelectrics. 

PACS numbers: 77.22.Ej, 77.80.-e, 77.80.Bh, 77.84.-s 
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Phase transformation phenomena are an important issue in sohd state physics. In un- 
derstanding these phenomena, both phenomenological and microscopical approaches have 
been playing comparable roles as one readily finds from the developments of the magnetism, 
superconductivity, and ferroelectricity. In particular, considerable attention is paid to the 
phenomenological Landau-Devonshire theory of ferroelectrics in the modern epoch where 
more and more advanced techniques of first-principle calculations of ferroelectrics become 
available. In the phenomenology of ferroelectrics, the issue recently attracting attention was 
the role of high-order anharmonic polarization terms in the Landau free-energy expansion 
In these works, it was demonstrated that the terms containing eighth (or higher) 
power of polarization, neglected in the traditional phenomenological framework, may be 
vitally important for the description the phase diagrams of perovskite ferroelectrics. The 
physics behind this effect is the need of high-order terms for adequate description of the sym- 
metry of the problem. Another demarche in the field has been recently undertaken by Li et 
al. , who have revisited the phenomenological theory of classical ferroelectric BaTiOs j^. 
These authors suggested incorporation of the eighth power terms in the Landau expansion of 
BaTiOs. In contrast to Refs. , here no symmetry arguments have been involved. The 
message by Li et al. was that the Landau expansion containing eighth power terms (with 
temperature independent anharmonic coefficients) is as efficient in the description of many 
properties of BaTiOs as the traditional sixth order expansion with the temperature depen- 
dent anharmonic terms jj,!^. This conjecture is of the interest since the strong temperature 
dependence of the anharmonic coefficients in the BaTiOs Landau expansion is in confiict 
with the displacive nature of ferroelectricity in the material, at least under the common 
assumption that the critical fiuctuation are weak in this system ~\ . Interestingly, in their 
analysis of the Landau potential using first principle calculations, Iniguez et al. obtained 
that the sixth order expansion accounts for the main features of BaTiOs phase diagram, but 
that all coefficients in the expansion have nontrivial temperature dependence j^. 

Together, these results pose two questions. First, are higher than sixth power terms in 
the BaTiOs Landau expansion needed for the correct physical description and are not just 
a matter of convenience? Second, does the incorporation of such terms enable elimination 
of the strong temperature dependence of the anharmonic coefficients in the BaTiOs Landau 
expansion? The present paper addresses these questions using the authors' experimental 
data and those available in the literature to unambiguously obtain a positive answer to the 
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first question and a negative to the second. 

To demonstrate the crucial role of the eighth-power term in the Landau expansion we will 
address the dielectric non-linearity of BaTiOa with respect to the electric field E applied 
along the [001] axis, in the vicinity of the cubic/tetragonal phase transition temperature. 
In this case only [001] component P of the polarization is involved, so that the Landau 
expansion of the Gibbs potential with the eighth-power term reads 

AG = ^-aP' + \(3P' + + l^P' (1) 

Accordingly we have for the equation of state and c-axis permittivity e^. 

E = aP + (3P^ + 7P^ + (2) 

e-^ = a + 3(3P^ + 5-fP* + 76P^ (3) 

These equations provide description of the non-linear dielectric response of BaTiOs in terms 
of the eighth-power expansion. One can compare this description with that in terms of the 
sixth-power expansion 

E = aP + (3'P^ + (4) 

e-^ =a + 3/3'p2 + syP' (5) 

which has been traditionally used in the field. These descriptions are clearly different. 
Formally, as one can readily check, the situation described by the expansion Eq. Q with 
the polarization-independent coefficients (3, 7, and 6 corresponds to polarization-dependent 
coefficients /?' and 7' in Eq. (jlj): 

f3' = (3- 6P\ 7' = 7 + 26P^ (6) 

In general, the consideration of the polarization-dependent anharmonic coefficients in the 
Landau expansion has no physical sense. However, if a not too wide interval of polarization 
variation is experimentally addressed with a finite measurement accuracy, coefficients (3' 
and 7' defined by Eq. (jH)) can be considered as polarization-independent when fitting the 
dielectric non-linearity. In other words, serious problems may arise in situations when a large 
variation of polarization is involved, e.g. the "jump" of the polarization on crossing the first 
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FIG. 1: (Color online) Schematic of the phase transition induced by an electrical field applied along 
[001] axis above Tc for BaTiOs single crystal. In paraelectric phase the polarization varies with 
field below Pq, while in ferroelectric phase the polarization is dominated by Pq = Pq + Ps- Curves 
(a) and (b) schematize the typical profiles of the Gibbs potential in paraelectric and ferroelectric 
phase, respectively. 



order phase transition. To verify this point we consider the case just above Tc, where the 
ferroelectric phase can be induced by an electrical field. As is schematized by the E ^ 
segment of the "double hysteresis loop" in Fig. Q the polarization "jumps" at the critical 
field Ec, whereas below or above this value the polarization varies slowly. At electric fields 
corresponding to the paraelectric phase, the contributions of the polarization dependent 
terms can be estimated from Eq. © to be not more than 5Pc^ and 25 Pc^, where Pc is the 
maximal values of the polarization in this phase. Using the 6 value of 2.9x10^^ Vm^^C"^ 
and the approximate value of Pq in our measurements (~0.04 C/m^), the magnitudes 
of the correcting terms in Eq.® are around 8x10^ Vm^C~^ and 8x10^ Vm^C^^ for f3 and 
7, respectively. In comparison with the typical values of (3 (-6~9xl0® Vm^C~'^) and 7 
(1^2x10^" Vm^C~^), these corrections are negligible to within at least of a few percent. 

At electric fields corresponding to the ferroelectric phase, the polarization is much larger 
than in the paraelectric phase. Provided that the interval of the polarization variation 
addressed is much smaller than the minimal polarization in the induced ferroelectric phase, 
Pq, when evaluating (3' and 7', one can take into account only the contribution of Pq in the 
correcting terms in Eq. ©, leading to the polarization independent (3' and 7' 

(3'^(3-6{Po)\^'^^ + 25{Pof. (7) 
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Using the typical value of Pq (~0.16 C/m^), the second terms in these relations can be esti- 
mated as 2x10*^ Vm^C~^ and 1.5x10^° Vm^C~^ for (3 and 7, respectively. These corrections 
are comparable with the typical values of the corresponding coefficients given above and 
cannot be neglected. 

Thus, from the above analysis, we conclude that if the eighth-power term in the Landau 
expansion exists having a value of the order of that estimated by Li et al. 0| , the anharmonic 
polarization coefficients (3' and 7' of BaTiOs estimated in terms of the sixth-power expansion 
should exhibit appreciable jumps on crossing the field induced phase transition. Such "phase 
sensitivity" of the nonlinear coefficients of the sixth order expansion enables the detection 
of the eighth-power term, suggesting a simple and direct way to assess the anharmonicity of 
BaTiOa. 

Our experiments unambiguously detect this phase sensitivity. We have performed pre- 
cise measurements of the dc field dependence of the dielectric permittivity of a BaTiOa 
single crystal slightly above Tc along [001] direction. An undoped BaTiOs single crystal 
(commercially available, GB group. Inc.) is cut and polished into a 2x2x0.3 mm'^ plate, 
with the main surface perpendicular to [001]. The two main sides are then entirely coated 
with platinum electrode to form a capacitor. A computer controlled setup, consisting of a 
temperature chamber, a high voltage power supply, and a precision LCR meter, was used 
to control or measure the temperature, the bias voltage, and dielectric properties (10 kHz), 
respectively. 

Figure El shows the experimental data measured at 135°C, which is around 4 °C higher 
than Tc- Little hysteresis effect related to the relaxation phenomena is found within either 
phase, indicating the high insulating quality of the crystal and the space charge free state 
of the crystal-electrode interface j^. The sixth order expansions, Eq. (jlj) and are used 
to fit the curves with the common value of a, and remarkably different (3 are obtained 
for either phase, as listed in Table HI As for 7, we cannot determine it only from the 
nonlinear dielectric behavior of the paraelectric phase. The position of Ec is used to obtain 
an estimate, which remarkably differs from that derived in the ferroelectric phase. The 
ferroelectric-phase coefficients have a serious problem in describing the dielectric properties 
of the paraelectric phase, and vice versa. As illustrated in the inset of Fig. |2l the eighth 
order expansion with a common set of coefficients enables a good fit in both phases at the 
same time. As-estimated values of /5, 7, and 5, agree well with the set suggested by Li et al. 
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FIG. 2: (Color online) Field dependence of the dielectric permittivity for BaTiOs single crystal 
at 135°C. The open squares represent the measured values, and the dotted line is to guide the 
eyes. Remarkably different (3 are obtained for either phase when fitting with Eq. (jlj) and 1(H)). 
as listed in Table ^ The dashed and solid lines are calculated with the coefficients derived from 
the paraelectric and ferroelectric data, respectively. The "ferroelectric" coefficients have serious 
problem in describing the dielectric properties of the paraelectric phase, and vice versa. The inset 
shows much better fit quality of an eighth order expansion for both phases. 



TABLE I: Nonlinear coefficients obtained from the field dependence of dielectric permittivity shown 
in Fig. m 



Order Phasc(s) -/3(Vm5C" 


-3) 7(Vm9C-5) <5(Vmi3c-'^) 


6 Paraelectric 8.6 X 10* 


1 X IQi" " 


6 Ferroelectric 18.7 X 10* 


8.2 X 10"' 


8 Both 8.6 X 10* 


0.8 X 10^0 3.4 X 10" 



"The value of 7 cannot be determined only by fitting the paraelectric branch. The present value is estimated 
by combining the dielectric nonlinearity and the position of the critical field Ec- 

a. 

Additional evidence for the eighth order term can be obtained using the dielectric per- 
mittivity values at Tc, where the phase transition is induced by temperature in the absence 
of external field. Taking into account that at Tc, E = Q and AG = 0, Eqs.Q and (jH)) lead 
to a simple relation for the variation of the dielectric permittivity: R = = 4 + 
where e^c and e^T stand for lattice permittivity along [001] at Tc in the cubic and tetrag- 
onal phases, respectively. One sees that, if 6 were zero or small enough, the ratio R would 
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TABLE II: Dielectric permittivity of BaTiOs single crystal at the phase transition temperature, 
Tc- 



Authors, Ref. ec,c ^c,T R 

Merz, [12] 14,200 2,300 6.2 

Drougard, [13] 16,000 2,200 7.3 

Meyerhofer, 16,500 5,000 3.3 

Johnson, [15j 9,130 1,500 6.1 

This work 11,400 2,200 5.2 

be close to 4 In Table |nj our data for R and those available in the literature are 

summarized. Apparently most of the values are remarkably larger than 4. Considering the 
possible experimental inaccuracy in determining the lattice permittivity in the ferroelectric 
phase due to not fully eliminated domain contribution, the real value of R is expected to 
be even larger than the presented values. The deviation of R from 4 provides additional 
evidence for the important role of the eighth order term. 

To elucidate more clearly the contributions of the various terms to the inverse dielectric 
permittivity, these are plotted as functions of polarization in Fig. IHl following Li et al. 's 
constants. It is seen that, in the low polarization regime, typically in the paraelectric phase, 
the quartic term dominates the whole nonlinear contribution and the dielectric nonlinearity 
can therefore be well fitted even by fourth order expansion. In the high polarization regime, 
typically in the ferroelectric phase, the contribution of the eighth order term becomes com- 
parable to or even larger than that of the sixth order term. 

Thus, we have shown that the eighth-power term is vitally important in the thermody- 
namics of BaTiOs, supporting the first conjecture by Li et al. . However, their second 
conjecture, that the eighth-power term can be introduced to avoid the use of temperature 
dependent anharmonic coefficients in the BaTiOs Landau expansion, remains unchecked. 
The direct way to do that is to analyze the dielectric non-linearity in the paraelectric phase 
in terms of the full expansion. As mentioned above, in the paraelectric phase, the dielec- 
tric nonlinearity is predominately controlled by the quartic coefficient (3. For this reason, 
the experimental results on a substantial temperature dependence of (3' of the sixth-power 
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expansion in the paraelectric phase [10|, llJ, lla ll7| strongly suggests that j3 coefficient from 
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FIG. 3: (Color online) Contributions of the various terms to the inverse dielectric permittivity, 
following Li et al. 's constants. In the low polarization regime, typically in paraelectric phase, the 
nonlinear contribution is dominated by the quartic term. In the high polarization region, typically 
in ferroelectric phase, the contribution of the eighth order term is comparable to or even larger 
than that of the sixth order term. 



the full expansion is also temperature dependent. However the aforementioned experimental 
data are rather scattered. To clarify this point, we measured the dielectric permittivity as 
a function of the field at a series of temperature points to find a reversal of the sign of 
the derivative at some 35 K above Tc (Fig. |^, implying a change of sign of /3 1^. Such 
behavior is clearly impossible on the assumption of temperature independent f3, even in 
terms of the eighth-power expansion. From this, we can definitely conclude that the ther- 
modynamic framework with a temperature independent /3 coefficient does not provide an 
adequate thermodynamic description of BaTiOa even in terms of the eighth-power Landau 
expansion. Thus, our results do not support the second conjecture by Li et al. . 

In conclusion, analyzing the non-linear and linear dielectric data on BaTiOs single crys- 
tal at temperatures close to ferroelectric/paraelectric phase transition we have demonstrated 
that the eighth-power anharmonic polarization term plays an essential role in the thermo- 
dynamics of the material. Our result justifies the corresponding conjecture by Li et al. 
3|. At the same time, based on the result on the temperature dependence of the dielectric 
non-linearity in the paraelectric phase, we have demonstrated that even the eighth-power 
Landau expansion should contain temperature dependent anharmonic polarization terms in 
order to provide an adequate thermodynamic description of the material. All these findings 
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FIG. 4: (Color online) Field dependence of dielectric permittivity at three representative temper- 
atures, evidencing the reversal of the sign of /3 at around 165°C 



attests to anomalously strong polarization anharmonicity of BaTiOs, which is unexpected 
for classical displacive ferroelectrics. 
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